Giant amplification of modes in parity-time symmetric waveguides 
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Abstract 

The combination of the interference with the amplification of modes in a waveguide with gain and losses can result in a 
giant amplification of the propagating beam, which propagates without distortion of its average amplitude. An increase 
of the gain-loss gradient by only a few times results in a magnification of the beam by a several orders of magnitude. 
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|1 Introduction 

0- . . 

• j— I • Amplification of guided optical waves is necessary to 
O j compensate for unavoidable losses along the propagations 
O -distance. The amplification is also important in the initia- 
te "tion process for the pulse generation requiring high intensi- 

. '^ es withj possibly, very low intensity inputs. Additionally, 
jthe concept of the gain-guidance, where the optical fibers 
^ Xand planar waveguides with gain allow for the efficient 
^ ^propagation of low power modes, was suggested recently [l[ 

1— i.as an alternative to the conventional index-guiding struc- 
1 jtures. The phenomenon was observed in the Nd-doped op- 

.tical fibers with low refractive index profiles [2j and served 
| as a basic element for the single-mode optical laser 
Q\ .Later on, the idea ofgain-guidance was extended also to 
the nonlinear media 0, Q ■ 

The gain element causes either the unbounded growth 
of the signal intensity or, if nonlinear dissipation becomes 
dominant, the convergence to some stationary or pulsating 
^ solution (i.e. to an attractor) 0]. In all cases the behavior 
• • -of the system is fully determined by the values of the sys- 
e ^ jtem parameters and does not depend on the characteristics 
^ .of the input signal (provided that the input signal belongs 
$_i "to the attractor basin). Therefore, in the linear systems 
^ one usually considers distortions of modes due to the gain 
(implemented, say, by doping), rather than propagation of 
such modes for long distances. 

On the other hand, a large wave intensity in some spa- 
tial domains can be achieved by the constructive interfer- 
ence of modes. In such amplification scenario, the increase 
of the intensity of the superimposed modes is determined 
solely by the characteristics of the input beam, and as such 
is limited by an a priori given value. 

In this Letter we propose to exploit the interference of 
the gain-guided modes which, on the one hand, allows for 
a flexible control of the output intensity by variation of 
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the input amplitude, and, on the other hand, benefit from 
the amplifying properties of the medium, thus resulting in 
a giant amplification of input beams over relatively short 
propagation distances. Moreover, by a proper design of 
the dissipation and gain profiles it is possible to achieve 
the guidance of modes for an arbitrary distance, without 
overcoming an a priori given intensity limit, and to realize 
cither single or multichannel guidance (selected by chang- 
ing the input beam) inside the same waveguide without 
changing the waveguide properties. Below, we focus on 
the linear case, when the effects of nonlinearity can be 
neglected. 

Our idea is based on the use of the special design for 
the gain and loss distribution to assure the parity-time 
(or VT) symmetry of the system. Suggested in the quan- 
tum mechanical context [7|, this concept implies spatial 
(V) and temporal (T) symmetry of the system, which in 
optics this can be implemented by means of even distribu- 
tion of the real part of the refractive index of the guid- 
ing medium combined with the anti-symmetrically dis- 
tributed gain and losses Particular relevance of the 
■PT-symmctric media for optical applications is justified 
by the property that in a definite range of parameters of 
such media all linear modes have purely real propagation 
constants (which corresponds to the pure real spectrum of 
the quantum- mechanical Hamiltonian Q). In other words, 
this means that all linear modes propagate without infi- 
nite growth or decay, thus allowing for the mode guidance 
in media with gain and losses. It is relevant to notice 
here that the VT symmetry may be not sufficient for ex- 
istence of purely real spectrum. The situation when a 
PT-symmctric system acquires (for some values of the pa- 
rameters) complex eigenvalues (propagation constants) is 
usually referred to as spontaneous PT-symmctry break- 
ing. In this Letter, however, we consider a situation when 
VT symmetry breaking does not occur. 
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The setup for the experimental implementation of PT- 
symmetric systems with gain and losses using the wave 
guiding structures was proposed in H, and later devel- 
oped in numerous studies (e.g., for the linear wave prop- 
agation see 0, [!(}). The first experimental studies of the 
optical "PT-symmetric structures were reported in [ll, 12 1. 
In particular, such phenomena as spontaneous PT sym- 
metry breaking and power oscillations violating left-right 
symmetry, have been experimentally observed. It was also 
suggested that appropriately designed combinations of the 
gain and losses in a dual core coupler can be used for the 
amplification of light signals [l3j and for the noise filter- 
ing 0. 



The left and right eigenvectors are mutually orthogonal 
and normalized: 

(0 n ,^m) = J </>n(r)0 m (r)dr = 6 numi 5 n2 ^ m2 , (5) 

where 5 n ^ n is the Kronecker delta. Therefore, for a given 
input <7o(r) = ?(0, r), the field evolution can be found in 
the form 

<z(z,r) = V c n e^0 n (r), c n = «,g ). (6) 

* — 

However, the right eigenmodes </> n (r) [as well as the left 
ones <^ n (r)] are not orthogonal in the usual sense. Their 
scalar product reads 



2. Some properties of "PT-symmetric potentials 

We consider propagation of a paraxial beam governed by 
the dimcnsionless field q obeying the Schrodinger equation 
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Hq, 



H 



-V 2 



V(r)+iW(r), 



(1) 



where r = (x, y), V = (d x ,d y ), r = |r|, V(r) describes 
radially symmetric modulation of the waveguide refractive 
index, while W(r) describes a spatial distribution of the 
gain and losses. We focus on the case 



V(-r) = V(r) and W(r) = -W(-r), 



(2) 



which is a necessary condition for the operator % to be 
PT symmetric. In particular, we choose |15| 



V{r)=r 2 , W(r) = 2g x x + 2g 2 y 



!g ■ r 



(3) 



where g = (31,32) = characterizes the gradient of 

the gain-loss term having the amplitude 3: g 2 = gj+g 2 - It 
is relevant to notice that the beam amplification reported 
below, although less pronounced because of the interfer- 
ence of a smaller number of modes, can also be observed 
in the one-dimensional setting. The two-dimensional set- 
ting, besides giving a considerably stronger amplification, 
offers some additional effects, such as the beam splitting, 
also considered below. 

For the stationary modes q(z, r) = e^ nZ </> n (r), we obtain 
the eigenvalue problem H<fi n = —/3 n <fin, where n = (n\, n 2 ) 
is an ordered pair of nomiegative integers identifying the 
guided modes. This problem has a pure real spectrum (lij : 
/3 n = — (2ni + 2n 2 + 2 + g 2 ) with the respective eigenfunc- 
tions given by 



6„(r) = 



Hni (x + igi)H ri2 (y + ig 2 ) 
% /2™i+"2ni!n 2 !7r 

-il(x+i gi ) 2 + (y+ig2) 2 ] 



Xf 



(4) 



where H n (^) are the Hermite polynomials. Since all /3 r 
are real, all the liner eigenmodes propagate undistorted. 

In addition, <fi n — </>* (the asterisk stands for the com 
plex conjugation) is an eigenvector of the Hcrmitian con 
jugated operator 'W: W^iftn = —fl n <j>n- The two sets, <£ r 



and (fi n , constitute the left and right complete bases [16 1. 



(<fin,<fim) = j (fi* n (r)(fi m (r)dr = D nimi (gi)D n2m2 (g 2 ), (7) 

which means that the optical energy distribution between 
the modes changes during the propagation. The energy 
transfer between the n-th and m-th eigenmodes is de- 
scribed by the functions 

D nm ( 9j ) = (±i flj .) n+m - 2 "24e^ (^) ±2 L^(-2g 2 ), (8) 



here (£) are the generalized Laguerrc polynomials, 



(0), 



L 
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signs " 

n, respectively. 



and " 



k = \n — m\, /i = min(n, m), and the 
" stand for the cases m > n and m < 
Since the Laguerre polynomials do not 



possess negative roots [l7j], the coefficients D nm (gj) 7^ 0, 
unless 3j = 0. Since the spectrum /3 n is equidistant, we 
conclude from © that the solution q(z, r) as well as the 
energy flow U(z) = J |g| 2 <ir are 7r-periodic with respect 
to z. We emphasize that though the equidistant spectrum 
simplifies considerably the analysis, it is not necessary for a 
giant amplification induced by the interference of the gain- 
guided modes, which can be observed also for other profiles 
of the refractive index. Notice also that energy oscillations 
are known to be a typical feature of 'PT-symmetric systems 
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We intend to explore the fact that, the modes which in 
the "conservative" (i.e. index guiding) waveguides have 
zero intensity at the center of the waveguide, for g 2 ^ 
carry nonzero power, and moreover, for a sufficiently large 
g 2 the respective modes of the PT-symmetric operator H 
may acquire the global maximum precisely at at r = 0. 
Moreover, we show that in the case of a PT-symmetric 
profile with a sufficiently large g 2 , the intensity maximum 
is shifted to the mode with a high quantum number n, i.e. 
with a large value of the propagation constant. This fea- 
ture leads to a dramatically different interference pattern 
for such modes, allowing for a constructive interference of 
the maxima of many different modes precisely at r = 0, 
thus strongly enhancing the intensity of the beam. At the 
same time, a destructive interference of the modes outside 
the beam center may fully eliminate all of the secondary in- 
terference fringes observable in the conventional patterns. 
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3. Giant amplification 

Let us consider the input beam of a Gaussian shape, 
which is the lowest eigenmode of a conservative waveguide 
with the same refractive index [i.e. of (HJ) — with g = 0]: 
90 W = v/t/oAe-^ 2 / 2 , where U = U(0) = J \q \ 2 dr is 
the input beam energy flow. A remarkable fact is that the 
respective solution q(z, r) can be found in an explicit form: 



^0 e -i(2+£/ 2 )2-2ig-rsin 2 z-ig 2 [sin(4z)/4-sin(2z)] 



xe 2g 2 sin 2 z g - l[( x - gi sin{2z)) 2 + (y-g 2 sin(2z)) 2 ]_ 



(9) 



Moreover, formula © is a particular case of a more general 
family of explicit solutions to ((H)-© given by [l8j |: 



q(z r) = y/Uo e i ' 3 n 2 +'' a ' rcos ( 2z ) _ia ' asin ( 4z )/ 4 

2a.gsi„ 2 z-g 2 /2 ^ (x - Xj (z) , y ~ X 2 (z)) 



xe 



^/L ni (-2gl)L n2 (-2g 



(10) 



where a = (01,02) is an arbitrary vector, which con- 
trols the angle of the incidence of the input beam, and 
X\ y 2{z) = ai,2 sin(22:) determine the position of the center 
of the mode. Indeed, © is obtained from (|10p simply by 
setting a = g and n = (0,0). The energy of the beam de- 
scribed by (TTU)) reads U(z) = C/oe 4a gsln z . In accordance 
with the expansion ©, U(z) is a 7r-periodic function. If 
the product a ■ g is positive, i.e. the beam is tilted to 
penetrate deeper in the domain with the gain, then U (z) 
approaches its maxima U max = U$e s at Zk = \ + 
and there is an exponential amplification of the input am- 
plitude achieved due to a combined effect of the gain-loss 
gradient g and a properly chosen direction of the input 
beam, i.e. a. Negative a • g, on the other hand, results in 
a periodic attenuation of the beam. 

Returning to expression ([9]) we find that the ampli- 
tude of the solution has the Gaussian shape with the cen- 
ter oscillating along the direction indicated by g. For 
< z < tt/2 the center of the beam is situated in the 
domain of gain, and the energy U(z) grows. The max- 
imal amplification is achieved in the point z = tt/2, i.e. 
when the center of the beam returns to the origin (where 
the gain is compensated by losses). For tt/2 < z < tt the 
beam resides in the lossy domain and loses the energy. At 
z = tt it returns to the origin and its energy approaches 
the lowest value Uq. Then the process is repeated (this 
recursive behavior is a consequence of the VT symmetry 
and resembles dynamics in a Hamiltonian system). 

To illustrate the above findings we focus, on the case 
of g = (51,0). At first, we consider the particular case 
given by The beam energy flow U(z) and the beam 
profiles at the maximal amplification point are shown in 
Fig. [U (a)-(c). Figure Q] reveals several important features. 
First, even a weak increase of the gradient results in a gi- 
ant amplification of the beam intensity as it is seen from 
panels (b) and (c) of Fig. [TJ Whereas these two panels are 
obtained for the same Gaussian input beam, an increase 
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Figure 1: (a) Energy flow vs propagation distance and (b) pro- 
file \q(ir/2, r)| 2 for the solution 10 with gradient g = (1/2,0). In 
the panel (b) the gain landscape W(r) for gi > and 32 = 
is schematically shown, (c) Profile \q(n/2, r)| 2 for the solution (0 
with g = (3/2,0). (d) Profile |g(7r/2,r)| 2 of the solution JTD) with 
n = (0, 1) and gradient g = a = (3/2, 0). In all panels Uo = 1. 



Table 1: Absolute values of the coefficients C( nii g) for Uo = 1, 32 
and different g\ [for 712 7^ all Cf„ lin2 j = 0]. 
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2 


3 


4 


9i 


= 1/2 


1.06 


0.37 


0.09 


0.02 


0.003 


9i 


= 3/2 


1.76 


1.86 


1.40 


0.85 


0.45 



of the gradient amplitude g by three times [from g = 1/2 
in panel (b) to g = 3/2 in panel (c)] results in magnifica- 
tion of the peak intensity by about 3 x 10 3 times [the peak 
intensity is about 0.8 in panel (b) against approximately 
2.5 x 10 3 in panel (c)]. Second, both in panels (b) and (c) 
the maximal beam intensity is achieved right at the center 
of the waveguide, although many eigenmodes are excited 
and, moreover, the maximal portion of the energy is con- 
centrated in a higher mode. To illustrate the last state- 
ment, in Table [1] we display the coefficients of the expan- 
sion ^ of the input Gaussian beam qo(r) = y/Uo/ne" 7 "^ 2 
over the eigenmodes (f) n : 



c n = v / C/o/(2»i+™ 2 n 1 !n 2 !) {i 9l )^ (* ff2 )" 2 e 9 <\ (11) 

From Table Q] one can see, that for g = (3/2,0) the most 
excited mode, i.e. the one having the largest coefficient 
c n , is the mode corresponding to n\ — 1 (rather than the 
mode with n\ = 0). With further increase of the gradient 
gi the maximal amplitude is shifted towards higher modes. 

The obtained results do not exhaust all possibilities of 
our setup. In particular, for the input beam with the 
shape of the first Gauss- Hcrmite function, the solution (fTOl) 
yields giant amplification of the split beam, as shown in 
Fig. [U (d). Even more sophisticated profiles of the ampli- 
fied beam can be obtained by variation of the input beam, 
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i.e. by changing n, or by using the incidence angle "ro- 
tated" with respect to the gain-loss gradient (i.e. using a 
different relation between g and a). 

4. Conclusion 

To conclude, we have proposed a physical setup where a 
giant but controllable, i.e. limited by the amplitude, am- 
plification of a guided beam can be achieved during the 
propagation. Such amplification occurs due to the inter- 
ference of the multiple gain guided modes. The latter is 
possible due to the reality of the propagation constants of 
all modes, on the one hand, and on the other hand due 
to the energy transfer among the modes which results in 
exciting many modes, practically by arbitrary input beam. 
The phenomenon is determined by the parameters of both 
the system and the input beam, thus allowing for an effi- 
cient managing over the output beam amplitude and, to 
some degree, over the output beam shape. This last prop- 
erty makes the configurations particularly relevant for the 
practical purposes whenever it is necessary not only to en- 
hance of an input beam until some given amplitude but 
also change this limit in the course of the experiment by 
managing only the input field. We also notice that further 
studies, including different 'PT-symmctric profiles allow- 
ing for the spontaneous symmetry breaking or nonlinear- 
ity, which may become relevant at high intensities of the 
phenomenon are of interest. 
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